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ABSTRACT 

Nonlinear diffusion equations of spectral transfer are systematically derived 
for anisotropic magnetohydrodynamics in the regime of wave turbulence. The 
background of the analysis is the asymptotic Alfven wave turbulence equations 
from which a differential limit is taken. The result is a universal diffusion-type 
equation in k-space which describes in a simple way and without free parameter 
the energy transport perpendicular to the external magnetic field B for trans- 
verse and parallel fluctuations. These equations are compatible with both the 
thermodynamic equilibrium and the finite flux spectra derived by Galtier et al. 
(2000); it improves therefore the model built heuristically by Litwick & Goldreich 
(2003) for which only the second solution was recovered. This new system offers 
a powerful description of a wide class of astrophysical plasmas with non-zero 
cross-helicity. 

Subject headings: MHD - solar corona - turbulence 

1. Introduction 

The observations of astrophysical plasmas by various spacecrafts have added substan- 
tially to our knowledge of magnetohydrodynamic (MHD) turbulence. Among the different 
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media widel y analyzed like the interstellar medium (jScalo and Ermegreenll2004l ) or the Sun's 
atmosphere (jChae et al.lll998l ). the solar wind is certainly the most interesting plasma since 
direct measurements are possible. This unique situation in astrophysics allows us to probe 
deeply the nature of the fluctuations and to investigate for e xample the origin of anisotropy 



( iMatthaeus et al. 



1996 



Alexakis et al.l 120071 ; iPodestal 120091 ) , to e valuate the mean energy 



dissipation rate (IMacBride et al.l 120081 ) . to detect intermittency (ISalem et al.l 120091 ). or to 



analyze the transition to the regime of dispersive turbulence characterized by a steepening 
of the magnetic field fluctuations spectrum with a power law index going f rom —5/3, at 



frequencies lower t han 1Hz, to indices lying around — 2.5 at higher frequencies (lGaltierll2006 
Smith et al.lbood : iGaltierlhoOsl : ISahraoui et aDl2009h . 



The low solar corona provides a second interesting example where it is believed that 
MHD turbulence plays a central role in the dynamics and the small scale heating. For ex- 
ample, in active regi on loops spectrom eter analyses revealed non-thermal velocities reaching 
sometimes 50 km/s (jChae et al.lll998l ); this line broadening is generally interpreted as un- 
resolved turbulent motions with length scales smaller than the diameter of coronal loops 
which is about one arcsec and timescales shorter than the exposure time of the order of few 
seconds. Turbulence is evoked in the sol ar coronal heating problem since it offers a natural 
proce ss to produce small scale heating (IHeyvaerts and Priestl Il992l ; iGaltierl Il999l ; ICranmer 
2010l ). Weak MHD turbulence is now proposed has a possible regime for some coronal loops 
since a very small ratio is expe cted between the fluctuating magnetic field and the axial com- 
ponent (IRappazzo et al.l 120071 ). Inspired by the obs ervations and by re cent direct numerical 
simulations of three-dimensional MHD tu rbulence (iBigot et al.ll2008al ). an analytical model 
of coronal structures has been proposed ( Bigot et al. 2008bl ) where the heating is seen as 
the end product of a wave turbulent cascade. Surprisingly, the heating rate found is non 
negligible and may explain the observational predictions. 

A third example where MHD turbulence seems to be fundamental is given by the upper 
solar corona which makes a connection between the lower corona and the stationary solar 
wind. Observations reveal that the heating in this region affects preferentially the ions in the 
direction perpendicular to the mean magnetic field. The el ectrons are much cooler than the 
ions, with temperatures generally less than or close to 10 6 K (IDavid et al.lll998l ). Additionally, 
the heavy ions become hotter than the protons within a solar radius of the coronal base. 
Ion cyclotron waves could be the agent which heats the coronal ions and accelerates the fast 
wind. Naturally the question of the origin of these high frequency waves arises. Among 
different scenarios, turbulence appears to be a natural and efficient mechanism to produce 
ion cyclotron waves. In this case, the Alfven waves launched at low altitude with frequencies 
in the MHD range, would develop a turbulent cascade to finally degenerate and produce ion 
cyclotron waves at much higher frequencies. In that context, the wave turbulence regime 
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was considered in the weakly compressible MHD case at low- /3 plasmas (where /3 is the 
ratio between the thermal and magnetic pressure) in order to analy ze the nonlinear three- 
wave interaction transfer to high frequency waves (jChandranl 120051 ) . The wave turbulence 
calculation shows - in absence of slow magnetosonic waves - that MHD turbulence is a 
promising explanation for the anisotropic ion heating. 

MHD turbulence modeling is the main tool to investigate the situations previously 
discussed. Although it cannot be denied that numerical res ources have been significantly 
improved during the last decades (jMininni and Pouquetl 120071 ). direct numerical simulations 
of MHD equations are still limited for describing highly turbulent media. For that rea- 
son, shel l cascade models are c urrently often used to investigate the small scale coronal 
heating (jBuchlin and Vellil 120071 ) and its impact in terms of spectroscopic emission lines. 
Transport equations are als o used for example in the context of solar wind acceleration in 
the extended solar corona ( jCranmer and van Ballegooijenl 120031 ) . The ad hoc model is an 
advection-diffusion equation fo r the evoluti on of the energy spectrum whose inspiration is 
found in the original paper by iLeithl ( 119671 ). It is also a cascade model where the locality 
of the nonlinear interactions is assumed but where the dynamics is given by a second-order 
nonlinear partial differential equation whereas we have ordinary differential equations for 
shell models. 

In next Section, the origin of the Leith's model is discussed and in Section [3] the Alfven 
wave turbulence equations are reminded in the case of non-zero cross-helicity. In Section HI 
the differential limit is taken on the previous wave turbulence equations and the associated 
nonlinear diffusion equations for anisotropic MHD turbulence are systematically derived. 
Finally, a conclusion is developed in the last Section. 



Leith's model 



A theoretical understanding of the statistics of turbulence and the origin of the power 
law energy spectrum, generally postulated from dimensional considerations d la Kolmogorov, 
remains one of the outstanding problems in classical physics which continues to resist modern 
efforts at solution. The difficulty lies in the strong nonlinearity of the governing equations 
which leads to an unclosed hierarchy of equations. Faced with that situation different models 
have been developed lik e closure models i n Fourier space for hydrod ynamic and magnetohy- 
drodynamic turbulence (jKraichnanlll963l ; lOrszag and Kruskallll968l ). In the meant ime - and 
follow ing an approach often fruitful in radiation and neutron transport theory (IDavidson 
19581 ) - Leith introduce d the idea o f a diffusion approximation to inertial energy transfer 
in isotropic turbulence ( ILeithl Il967[ ). This new class of ad-hoc models describes the time 
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evolution of the spectral energy density, e(k), for originally an isotropic three-dimensional 
incompressible hydrodynamic turbulence, in terms of a partial differential equation by mak- 
ing a diffusion approximation to the energy transport process in the k-space representation. 
Ignoring forcing and dissipation the three-dimensional isotropic Navier-Stokes equations read 
in Fourier space 

The radial component of the energy flux vector is modeled as 

J> = -2>(*)^, (2) 

where D(k) is a diffusion coefficient that remains to be determined. It is straightforward to 
show dimensionally that the diffusion coefficient scales as 

_ k , „ 

(3) 

r 

where r is the typical transfer time of the Navier-Stokes equations which can be identified 
as the eddy turnover time r e ddy Therefore, we may evaluate this time as 

1 

T = T eddy ~ ; r — . (4) 



kv ek 3 

We remind that the total kinetic energy per mass is by definition J e(k)dk = J E(k)dk. 
After substitution of (jlj) into (J3]) it is possible to r ewrite (up to a factor) the model equation 



JTJ) for the omnidirectional spectrum ( iLeithl 119671 ) 

dE{k) 



dt 



(5) 



which is commonly named the Leith's equation. Beyond its relative simplicity, equation (151) 
exhibits several important properties like the preservation after time integration of a non 
negative spectral energy and the production of the Kolmogorov spectrum in the inertial 
range which corresponds to a finite energy flux solution. It is straightforward to prove that 
by imposing a constant energy flux in the inertial range, namely 

fc H/2£l/2 3 (£ /jfe 2) = constant . ( 6 ) 

If we look for power law solutions, E ~ k x , then the unique solution that emerges is x = —5/3. 
Note that this equation may als o exhibit an anomalous scaling du ring the non stationary 



phase with a steeper power law (jConnaughton and Nazarenkoll2004( ). 
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A generaliza tion of the Leith's model to three-dimens ional isotropic M HD turbulence 
was proposed by IZhou and Matthaeusl (119901 ) . (Note that llroshnikovl ( 119641 ) proposed the 
first such a model for MHD from which the —3/2 spectrum was derived). The main modifi- 
cation happens in the evaluation of the transfer time r for which a combination of the eddy 
turnover time r e ddy and the Alfven time ta is proposed. The phenomenological evaluation 
of the transfer time allows the recovering of either the Heisenberg-Kolmogorov (—5/3) or 
the Iroshnikov-Kraichnan (—3/2) spectrum when the ratio t py m ,, /t4 is respectively much less 
or m uch larger than one (jKolmogorovl 11941c iHeisenbergl Il948l ; llroshnikovl Il964t iKraichnan 
19651 ). The model was also adapted to the case of a non-zero cross-helicity for which a 
distinction was made between the Elsasser energies E + and E~. The generalization of the 
Leith's model to the more realistic situation of anisotropic MHD turbulence where an ex - 
ternal magnetic field B is imposed was proposed on ly recently (IMatthaeus et al.l 120091 ). 
As already announced by IZhou and Matthaeusl (Il990l ) the departure from the assumption 
of isotropic turbulence generates a difficult mathematical treatment since, in particular, a 
diffusion tensor is expected instead of a scalar. Another difficulty comes from the locality of 
the nonlinear interaction which is assumed in the isotropic case: when a mean magnetic field 
is imposed the situation is different since a reduction of nonlinear transfers occurs along B . 
In terms of triads, k = k + L, it means that one of the wa yevectors, say /t, is ma i nly or iented 
transverse to B . The sophisticated model proposed by IZhou and Matthaeusl (Il990[ ) is an 
attempt to describe such a nontrivial dynamics. 

The case of Alfven wave turbulence for which a re latively strong B n i s required is an im- 
portant limit for which a rigorous analysis is possible (IGaltier et al.ll2000l ). The wave kinetic 
equations derived are a set of coupled integro-d ifferential equations which are not obv ious to 
simulate numerically in the most general case (IGaltier et al.ll2000t iBigot et al.l 12008a ). This 
remark was a motivation for deriving a model made of two cou pled diffusion equations which 
describe Alfven wave turbulence with a non-zero cross-helicity (ILithwick and Goldreichll2003l ) . 
These model equations are abl e to recover the fin ite flux spectra which are exact solutions 
of the wave kinetic equations ( IGaltier et al.l 120001 ) . In the present paper, it is shown that 
a set of two coupled nonlinear diffusion equations may be derived systematically from the 
asymptotic equations of Alfven wave turbulence by taking a differential limit. An important 
difference is found between the non li near diffusion equations derived here and the model 
proposed by lLithwick and Goldreichl (120031 ) . The main physical problem is the inability for 
the model to reproduce the thermodynamic equilibrium solutions which are exact solutions 
of the wave turbulence equations. It is believed that the higher degree of accuracy of the new 
system offers a powerful description of a wide class of astrophysical plasmas with non-zero 
cross-helicity. 
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3. Asymptotic theory of Alfven wave turbulence 



The wav e turbulence theory for three-dimensional incompressible MHD was derived 
rigorously by iGaltier et al.l ( 120001 ) . It is a perturbative theory which necessitates heavy 
calculations which will not be reproduced here. We refer the reader to the original paper 
for a global explanation or to two satellite paper s where simplified approaches are adopted 



( iGaltier et al.l 120021 ; IGaltier and Chandran 



20061 ) . Since it is important to understand the 



wave turbulence equations from which our analysis will start, we shall remind below the 
main steps in their derivation. 



The inviscid incompressible three-dimensional MHD equations read 



d t z s - sB • Vz s 
V-z s 



• Vz s - VP* 



(7) 



0. 



where z s = v + sb are the Elsasser fields (s = ±), v is the fluid velocity, b is the magnetic 
field in velocity units, B is a uniform magnetic field (in velocity units, i.e. the Alfven 
speed) and P* is the total (thermal plus magnetic) pressure. We assume that the uniform 
magnetic field is relatively strong (Po ^ z s ) and that MHD turbulence is dominated by a 
wave dynamics for which the nonlinearities are weak. In such a limit, a small parameter e 
may be introduced formally to measure the strength of the nonlinearities. Then, we obtain 
for the jth-component 



(d t - sB Q d\\) 



Z 3 



—ez„ 



~ s d z s 



d P 



(9) 



where the Einstein's notation is used for the indices. Note that the parallel direction (||) 
corresponds to the direction along Bo. We shall Fourier transform such equations with the 
following definition for the Fourier transform of the Elsasser field components z|(x, t): 



z »(x, t) = J a*(k, t) e i ^ +SUJkt) dk , 



(10) 



where Uk = B$k\\ is the Alfven frequency. The quantity a|(k, t) is the wave amplitude in the 
interaction representation, hence the factor e lsuJkt . Then, the Fourier transform of equation 
© gives 



dtttjik) = -itk m P jn J J a m s (K)a s n (L)e lsAuJt S k:KL dKdL . 



Here, P jn is the projector on solenoidal vectors such that P, n (k) = 5j n — kjk n /k 2 ; 5k,KL = 
5(k — k — L) reflects the triadic interaction, and Au = ul — oju — u) K is the frequency mixing. 
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The appearance of an integration over wave vectors k and L is directly linked to the quadratic 
nonlinearity of equation (Q (as a result of Fourier transform of a correlation product). 

Equation f lTTj) is nothing else than the compact expression of the incompressible MHD 
equations when a strong uniform magnetic is present. It is the point of departure of the wave 
turbulence formalism which consists in writing equations for the long time behavior of second 
order moments. In such a statistical development, the time-scale separation, TA/r e dd y <C 1 
(with ta = 1/wfe and T e dd y = t/k±z s ), leads asymptotically to the destruction of some 
nonlinear terms - including the fo urth order cumulants - and only resonance terms survive 



(Galtier et al. 


2000; 


Galtier 


2009) 



equations. In such a statistical development, the following general definition for the total 
(shear- plus pseudo-Alfven wave) energy spectrum is used; 

(a*(k)a*(k')) = E s (k) 6(k + k')/k ± , (12) 

where () stands for an ensemble average and k± = |kj_|. In absence of magnetic helicity 
and in the case of an axially symmetric turbulence, the asymptotic equations simplify. For 
shear- Alfven wavea3, the energy spectrum is given by 

%»r(*±.*||)=»(*||)^(*±), (13) 

where g(k\\) is a function fixed by the initial conditions (ie. there is no energy transfer 
along the parallel direction). In the limit k± 3> k\\, the transverse part obeys the following 
nonlinear equation (the small parameter e is now included in the time variable and the limits 
of the integration are explicitly written) 

d t E s L (k ± ) = — / cos 2 0sinfl— E7 s (k l ) [k ± EUL ± ) - L ± E s ± (k ± )]dK ± dL ± , (14) 

B Jo Jo «L 

where is the angle between kj_ and Lj_, and 9 is the angle between k^ and K± with the 
perpendicular wave vectors satisfying the triangular relation kj_ = Lj_ + k± (see Figure [T]). 
Note that from the axisymmetric assumption, the azimuthal angle integration has already 
been performed, and we are only left with an integration over the absolute values of the two 
wave numbers, k± = \k±\ and L± = |Lj_|. In the same way, equations can be written for 
pseudo-Alfven waves which are passively advected by shear-Alfven waves, namely 

d t EUk ± ) = — / sintf— EI s (k ± ) [k ± ES(L ± ) - L ± EUk ± )] dn ± dL ± , (15) 

Bo Jo Jo K -L 



1 We recall that shear- Alfven and pseudo-Alfven waves are the two kinds of linear perturbations about 
the equilibrium, the latter being the incompressible limit of slow magnetosonic waves. 
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Fig. I.- 



Triadic interaction k = k + L and its projection in the plane perpendicular to B . 
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with by definition 

E^oik^k^^g^Elik^), (16) 
where g(k\\) is a function determined by the initial condition. 

Equations (fl~4l and ffT5]) are asymptotically exact. These master equations of Alfven 
wave turbulence describe the nonlinear evolution of MHD turbulence in the presence of a 
strong uniform magnetic field, with non-zero cross-helicity and zero magnetic helicityEl- In 
the limit k± 3> k\\ considered, equations (jHj) and (fT5|) describe respectively the dynamical 
evolution of transverse and parallel fluctuations. 



4. Differential limit for strongly local interactions 



We shall take a diff erential limit of equations (|14p and f JT5|) for strongly local interactions 
(jDyachenko et al.lll992l ). It is important to note that by taking this limit we shall extract a 
subset of the full set of interactions that is present in the Alfven wave turbulence equations 
(|14p - (|15p . In terms of triads, strongly local interactional means that we will only retain 
triangles which are approximately equilateral. (Note that the locality concerns only perpen- 
dicular wavevectors.) Therefore, the differential limit will lead to an approximate description 
of Alfven wave turbulence which is believed, however, sufficiently rich to reproduce the most 
important properties of the original system. The rigorous derivation will be presented only 
for shear-Alfven waves since the generalization to pseudo-Alfven waves is straightforward. 
Multiplying equation (j!4p by an unknown function f(k±) we obtain, after integration in k±, 



P + OO P + OO 

d t / E s ± (k ± )f{k ± )dk ± = / 

Jo JO 



where 



T(k_ L ,L_ L ,K±) 



71 



Jo 



kiL, 



cos 0sm< 



Note the use of the triangle relation (see Figure [1]) 

smu = sm0 — . 

K± 



T(k±, L±, K,±)f(k±)dk±d,K±dL_i 



EZ'(K ± )[k x E' ± (L ± ) - L ± E s ± (k ± )} 



(17) 



(19) 



2 We refer to the original paper (jGaltier et al.1120001 ) for a discussion about the domain of applicability in 
terms of wavevectors of the Alfven wave turbulence regime. 

3 Strongly nonlocal interacti ons is another intere sting limit from which we may deriv e turbulent viscositie s 
in the wave turbulence regime jBigot et al.ll2008bl) or in the strong turbulence regime |Pouauet et al.lll976l) . 
In the latter case, an EDQNM closure model was used and the main application was the dynamo problem 
whereas in the former case the application was the solar corona with coronal loops and coronal holes. 



It is convenient to introduce the following definition S = tt cos 2 <fi sin cf). Then, by changing 
the name of indices one can write 

dt / E s ± (k ± )f(k ± )dk ± = (20) 
Jo 

- / / [T(k±,L±,K±)f(k±) + T(L ± ,k±,K ± )f(L ± )}dk ± dK ± dL±. 

* Jo Jo Jo 

However, by symmetry we also have 

T{L ± ,k±,K±) = -T(k±,L±,K±) , (21) 

which gives 

1 



dt / E s ± (k ± )f(k ± )dk A 



oo ^+oo P + OO 



T(k ± , L ± , n ± )(f(k ± ) - f(L ± ))dk ± dK ± dL ± .(22) 



'0 JO 

For strongly local interactions, we have the following relations 

«± = k±(l + e K ) , |e K | < 1 , 

£± = fc ± (l + e L ) , |e L |<l, 
where e K and are two variables of small amplitude. Then, at first order we have 

df(k ± ) 



and also at first order 



f(L ± ) = f(k ± ) + e L k A 



S 



dkj 



(23) 
(24) 

(25) 



T(k ± ,k ± (l + e L ),k ± (l + e K )) = —E~ s (k ± )k ± e L 

-DO 



-^^r(^)fci ei ^- (jB1(A;±)Aj 



(26) 



Therefore, equation (|22|) may be reduced at first order as 



+oo 



d t / E s ± (k ± )f(k ± )dk A 





1 

-Bo 
1 

-Bo Jo 



+oo /•+£ r 



^Erik^klel-^iEKk^/kj 



dfjkj 
dki 



dk±de K deL 



C ± EZ s (k ± )kl-^(E s ± (k ± )/k ± 



dfjkj 
dk. 



dk i 



(27) 
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where 



C ± = J 6 J £ ^e 2 L de K de L} (28) 

with e a small positive number. The introduction of e is necessary to ensure the strong 
locality of the interactions and therefore the convergence of integral (|28|) . An integration by 
parts gives 



r+oo i iH 

dt / E> ± (k ± )f(k x )dk ± = — 
Jo B o Jo 



a 



dk 



J_ 



C ± EZ s (k ± )k 6 ± -^(E s ± (k ± )/k A 



f{k ± )dk ± . {29) 



Note that this operation implies for the function / some constrains of convergence at the 
boundaries. Since / is an arbitrary function one can write 

(30) 

which is the differential limit of the wave turbulence equation for shear- Alfven waves. It is 
useful to get an evaluation of C± since the constant in front of an equation always enters 
into account in the evaluation of the time scale dynamics. By noting that for strongly local 
interactions the angles of the triads are approximately tt/3, one obtains 

c -^£/->^=iS' (31) 

Clearly, the degree of locality will strongly modify the time scale. For example, for strictly 
local interactions e = and no evolution of the spectra is expected. This is consistent with 
the original equation (TH1) for which the right hand side is trivially zero if L± = k±. 



The same type of analysis for pseudo- Alfven waves gives 

a > E ^ = wA (tf™*: 1 1 • (32 > 

where C\\ = f_ e f_ e (ir/2)e 2 L sm<f)de K de L = 7ie 4 /\^3. 




5. Finite flux solutions and Komogorov constants 

Equations (150|) and (152"]) are the main results of the paper. They describe respectively the 
dynamical evolution of perpendicular and parallel fluctuations to the background magnetic 
field Bo- We see that in the differential limit of strongly local interactions the wave turbulence 
equations are much simpler. They still satisfy the finite flux solutions as we will see below 
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by looking at the power law solutions for shear-Alfven waves. First of all let us introduce 
the energy flux P±(k±) which is by definition 



We obtain 



d t E 8 ± (kj 



C, 



dPi(k A 



dk ± ' 
d fE s ± (k A 



We shall find the power law solutions by introducing E s ± = C s k± na into ( 135]) ; one gets 



Pl(k±) = -^C s C- s (n s - l)k 

-DO 



n s +n_ s +4 



with by symmetry 



C S C~ 



PIBo 



P7 S B 



Bn 



PfP7 



C ± {l-n s ) C ± (l-n_ s ) Cj_ V (l-0(l-n_ s ) ' 



which leads to 



Bn 



E±(k ± )E2(k ± ) = -± 



P+P7&7 4 



B n 



C± y 5 + 72+72 

4\/35 



PtPVk- 4 



1 



7re 4 y 5 + n + n_ 



-7,-4 



For balance turbulence, n+ = n. 



P±(A;j 



-2, and the finite flux solution is 

"Kkf 
■xfp^k- 2 



B 

3C ± X 
4P 



(33) 



(34) 



(35) 



Therefore, the finite flux solutions, P±{k±) = P± = constant, correspond to 

77 + +77_ = -4, (36) 



(37) 



(38) 



(39) 



13 



We see that the Kolmogorov constant depends on an arbitrary truncation of t he integration 



domai n in equation ( |3~T]) . We remind that the Kolmogorov constant found by iGaltier et al. 
( I2000l ) for balance turbulence was, Ck = 0.585-\/Po; therefore the constants coincide for 
e ~ 1.21. Although the previous value violates the assumption of smallness for e it could 
be taken a posteri ori to find a unique solution compatible with the exact derivation of 
Galtier et all (boooh . 



A similar analysis for pseudo-Alfven waves ( 1321) gives the energy flux 




(40) 



By introducing ES = C s kj_ ms we find the finite flux solutions m s + n„ 

m + + m_ = —4 , 

c + c- 



-4 and thus 



and 



E+(k ± )E7{k A 



5 + m + m_ ^p+p- 
B C ± V5 + n + n_ P{P{\ 



which reduces for balance turbulence to 

EMj 



5 + m + m_ ^p+p 



— n ± ' 



(41) 
(42) 

(43) 



'BnC, Pi 



3C, 



Bn 



Pi 



47rV3e 4 VP± 
0.214 ^ ^ fcT 2 



(44) 



Note that the Kolmogorov constant found by IGaltier et al.l (120001 ) for balance turbulence 
was, C' K = 0.0675a/Po; in this case the constants coincide for e ~ 1.78. 

Additionally, equations ( 13"Uj) and (1321 reproduce th e thermodynamic equilibrium solu- 
tions which correspond to zero flux ( IGaltier et al.ll2000l ). In this case, it is straightforward 
to show from (|30|) and (132)) that 

n s = m s = 1 . (45) 

As explain ed above, equations ([30]) and (|32|) are different from the model equations (73)-(74) 
derived in iLithwick and Goldreichl ( 120031 ) (where the notation are different; a comparison 
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is possible if one takes ke^ ~ E±~) which do not give the thermodynamic equilibrium 
solutions. The new system systematically derived here improved therefore the previous 
description while keeping the simplicity of a diffusion model. 

The nonlinear diffusion equations ( 130|) - (|32|) for non-zero cross-helicity may exhibit differ- 
ent power laws as finite flux solutions. Our knowledge of the initial system ( !T4])-ffT5l) imposes 



a pri ori that the power law indices satisfy the condition —3 < n s ,m s < —1 (IGaltier et al. 



20001 ) . However, if we look at the diffusion equations we do not find any other constrain 
than relations (I3"6"j) and (1411) which means that in principle the power law indices are not 
bounded. The simplicity of the diffusion equations allows us to write a simple relation for 
shear-Alfven energy fluxes, namely 

P+ Eld(El/k L )/dk L dln(Ef/k ± )/dk± 



Pi E+d(El/k ± )/dk ± d\n(Ei/k ± )/dk± 
In the stationary state, we obtain 

Pf n+ — 1 { n . — 



(46) 



P, ri- — 1 \ n 



(47) 



which gives P±/P± =1/2 for n + = — 1 and P±/P± = 2 for n + = —3. Note that the zero 
cross-helicity case corresponds to n + = n_ = —2 for which P_^ = PJ_*. In the general case 
which includes nonlocal inter actions, we remind that we have P±/P± = for n + = — 1 and 



P±/P± = +oo for n + = —3 (IGaltier et al.ll2000l ). Therefore, the differences found between 



both predictions (from the diffusion and the integro-differential equations) give an evaluation 
of the nonlocal contributions. 



6. Numerical illustrations 

In order to check if the constant flux solutions are attractive we have performed two 
numerical simulations of the nonlinear diffusion equations. Only the case of shear-Alfven 
waves (transverse fluctuations) has been considered. A linear viscous term is added in order 
to introduce a sink for the energy. In practice, the following equations are simulated 

dlEt(k±) = JL (SIM) ) - „ klEi{k±) , ,48) 

where v is the viscosity (a unit magnetic Prandtl number is chosen). This type of equations 
is favorable to the use of a logarithmic subdivision of the k± axis such that in our case 



k u = Tl 



(49) 
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where i is a positive integer. Such a discretization allows us to reach Reynolds numbers 
much greater than in direct numerical simulations. We take i max = 200 which corresponds 
to a ratio of about 10 6 between the largest and the smallest scales. In our simulations the 
viscosity is fixed to v = 5 x 10~ 5 . 

The first simulation corresponds to the zero cross-helicity case for which by definition 
E~\_ = E±. Large scale spectra centered around k± = fco are taken initially with the form 

E±(k ± )~k ± 3 eM-ki_ 2 /k 2 Q ), (50) 

where k = 5. Only the time evolution of E~*[ is given in Fig. [21 We see that the front 
of the energy spectrum propagates towards larger wavenumbers to reach eventually a kj 2 - 
stationary spectrum as predicted by the theory. We may note the acceleration of the front 
until the dissipation scale is reached since spectra are separated by a constant interval of 
time. In the second simulation we have fixed initially the (reduced) cross-helicity to 

p= ^r - 8 - (51) 

Like in the previous case the initial spectra are centered around k± = ko = 5 with more 
energy in spectrum E~\_ than in We keep the same form as (1501) . The result is shown in 
Fig. |3j We see that the compensated spectra fit well with the theoretical prediction fcj over 
several decades. In fact, spectra meet at relatively small k± and exhibit the same inertial 
range with the same fc_i_ _2 -spectrum over a wide range of scales. 



7. Discussion 



It would be releva nt to investigate if whether or nor the system recently derived by 
Matthaeus et al.l (120091 ) for strong (anisotropic) MHD turbulence is able to recover the 
present equations when the limit of wave turbulence is taken. It would also be interest- 
ing to analyze if an anomalous scaling is detected during the front propagation (which is 
not easy to find here). We remind that the wave kinetic equations (JHJ) exhibit a k±~ 7 ^ 3 - 
spectrum during the non-stationary phase (at zero cross-helicity) which is still not understood 
(IGaltier et al.ll2000f ). Ano malous scalings are weaker in di ffusion models of turbulence than 
in wave kinetic equations ( jConnaughton and Newell! 120101 ) . We plan to further investigate 
this point by using for example a higher order numerical scheme. We also plan to further 
compare numerically the nonlinear diffusion equations and the wave kinetic equations to 
dete rmine the domain of divergen ce between them or the influence of an external force (see 



e.g. 



Galtier and Nazarenkoll2008l ). 
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10° 10 1 10 2 10 3 10 4 10 5 10 6 



Fig. 2. — Time evolution of the energy spectrum E±(k±) for a cross-helicity p = 0. The 
stationary spectrum is well fitted by a power law in kj 2 . Forty spectra are shown with a 
constant interval of time. 
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10° 10 1 10 2 10 3 10 4 10 5 10 6 



Fig. 3. — Time evolution of the compensated energy spectra E^E±k± 4 for a cross-helicity 
p = 0.8. The stationary spectra satisfy the relation n + + n_ = —4. Forty spectra are shown 
with a constant interval of time. 
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The nonlinear diffusion equations for non-zero cross-helicity fl30l) - fl32|) is a simple and 
therefore useful system for describing a wide class of astrophysical plasmas. The solar corona 
with the myriad of magnetic loops which are characterized by a stro ng axial magnetic field 
is probably a good example of application of Alfven wave turbulence (IRappazzo et all 120071 ; 



Bigot et all l2008bl ) . This regime is also relevant for coronal holes where the solar wind 
is produced. For both examples, equations (!30|) - (l32|) could give a description of MHD 



turbulence at large scales since it seems inevitable that at smaller scales the strong turbulence 
regime overc omes. In this ca se a coupling with for example the advection-diffusion model 
proposed by IChandranl (120081 ) would be relevant. 



We acknowledge Institut universitaire de France for financial support. 
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